How to detect the fourth order cumulant of electrical noise 
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It is proposed to measure the current noise generated in a mesoscopic conductor by macroscopic 
quantum tunneling (MQT) in a current biased Josephson junction placed in parallel to the conductor. 
The theoretical description of this set-up takes into account the complete dynamics of detector and 
noise source. Explicit results are given for the specific case of current fluctuations in an oxide layer 
tunnel junction, and it is shown how the device allows to extract the fourth order cumulant of the 
noise from the MQT data for realistic experimental parameters. 

PACS numbers: 72.70.+m,05.40.-a,74.50.+r,73.23.-b 
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Within the last decade electrical noise has moved into 
the focus of research activities on electronic transport in 
nanostructures 0, since it provides information on mi- 
croscopic mechanisms of the transport not available from 
the voltage dependence of the average current. Lately, at- 
tention has turned from the noise auto-correlation func- 
tion to higher order cumulants of the current fluctuations 
characterizing non-Gaussian statistics 0,0 ■ While theo- 
retical attempts to predict these cumulants for a variety 
of devices are quite numerous experimental obser- 
vation is hard because of small signals, large bandwidth 
detection, and strict filtering demands. A first pioneering 
measurement by Reulet et al. of the third cumulant of 
the current noise from a tunnel junction has intensified 
efforts and several new proposals for experimental set- 
ups have been put forward very recently, some of which 
are based on Josephson junctions (JJ) as noise detec- 
tors. Lindell et al. Q employed a Coulomb blockaded JJ 
to demonstrate that the conductance of the junction in 
the Coulomb gap region is sensitive to the non-Gaussian 
character of noise applied to the junction. A modification 
of this set-up was suggested by Heikkila et al. @ to get 
specific information on the third cumulant of the noise. 
Another recent experiment has observed activated- 
over-the-barrier-jumps of a JJ biased by a noisy current. 
The data are consistent with resonant activation pro- 
duced by the second cumulant of the noise at the plasma 
frequency of the junction. For a measurement of the full 
distribution of current fluctuations Tobiska and Nazarov 
@ suggested to use an array of overdamped JJs acting 
as a threshold detector for rare current fluctuations trig- 
gering over-the-barrier jumps. Since in the overdamped 
limit retrapping spoils the built-up of a detectable volt- 
age, it was argued by Pekola that an experimentally 
more accessible detector would extract the noise charac- 
teristics from modifications of the macroscopic quantum 
tunneling (MQT) rate in an underdamped JJ. 

In all experimental set-ups to measure higher order 
cumulants realized and proposed so far, heating is one 
of the major experimental obstacles |Toj . Thus, experi- 
ments have primarily attempted to establish just the un- 
specified non-Gaussian nature of the noise or to measure 
the third cumulant (skewness). The latter one is particu- 
larly accessible since it can be discriminated from purely 
Gaussian noise due to its asymmetry, e.g. when inverting 
the current through the conductor. This is in contrast to 
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FIG. 1: Electrical circuit containing a mesoscopic conductor 
G in parallel to a JJ with capacitance Cj and coupling energy 
Ej biased by an external current 1^. The switching out of the 
zero voltage state of the JJ by MQT is detected as a voltage 
pulse V. 



the fourth order cumulant (sharpness), which on the one 
hand due to heating effects may be completely hidden be- 
hind the second and the third one, but on the other hand 
is required to gain an essentially complete characteriza- 
tion of the distribution of current fluctuations. In this 
Letter we propose and analyze a set-up, with the circuit 
diagram depicted in Fig. ^ which allows to detect the 
fourth order cumulant of the current noise generated by 
a nanoscale conductor. Since this conductor is placed in 
parallel to a current biased JJ in the zero voltage state, no 
heating occurs prior to the decay of this state by MQT. 
However, the MQT rate is modified in a specific way by 
the even higher order cumulants characterizing the non- 
Gaussian current fluctuations of the conductor. 

The complete statistics of current noise generated by a 
mesoscopic conductor can be gained from the generating 
functional 
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where I(t) is the current operator and T the time order- 
ing operator along the Kadanoff-Baym contour C. Time 
correlation functions of arbitrary order of the current are 
determined from functional derivatives of G[<j>], in partic- 
ular, the average current 

c 1 (t) = (i(t)) = ie ds G [4>]/dm^=o 
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and the current auto-correlation function 

c 2 (m') = (i(t)i(t')) = e 2 d?s G [<f>}/d<t>(t)d<t>(t')\ 



0=0 



Higher order functional derivatives give the cumulants 
related to non-Gaussian current fluctuations 
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We remark that the functional Sg [</>] carries the full fre- 
quency dependence of all current cumulants and not just 
their time averaged zero frequency values usually studied 
in the field of full counting statistics . 

By way of example let us consider an Ohmic resistor 
of resistance R in thermal equilibrium at inverse temper- 
ature (3. Then, the functional Sg[4>] = Sr[4>] takes the 
well-known form 

Sn[<f>] = \^Jl j c dt j c dt ' a (* - *0 ^ 



where 
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The quadratic form reflects the Gaussian nature of the 
current fluctuations in this case which implies that all 
cumulants except for C<z vanish. On the other hand, 
for a tunnel junction with many transmission chan- 
nels, where each channel has a small transmission co- 
efficient Tj leading to the dimensionless conductance 



g T = h/{ATTe 2 R T ) = 
ing resistance, one has 



where Rt is the tunncl- 



S T [4>] = -Ag T dt d£ a(t - t') shr 
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Here, the periodicity in <j) reflects the discreteness of the 
transferred charges associated with non-Gaussian current 
fluctuations. 

To gain information on the noise of the conductor, it 
may be placed in parallel to a current biased JJ as de- 
picted in the circuit diagram of Fig.^ For a bias current 
lb below the critical current I c , the J J is in its zero voltage 
state and the bias current flows as a supercurrent entirely 
through the J J branch of the circuit. Consequently, no 
heating occurs in the conductor and the total system can 
easily be kept at low temperatures, where the decay of 
the zero voltage state occurs through MQT. The rate of 
this process depends with exponential sensitivity on the 
current fluctuations of the conductor so that the JJ acts 
as a noise detector. 

The MQT rate T can be calculated in the standard 
way |12L If3| from the imaginary part of the free energy 
F, i.e., 

T = | lm{F} , 

where F = — m(Z) is related to the partition func- 
tion Z = Tr{e~P H }. In the path integral representation 
one has 



Z = J V[0]e 



S[9] 



which is a sum over all imaginary time paths with period 
h(3 of the phase difference 9 across the JJ weighted by 
the dimensionless action S[9] = S,jj[9] + Sg[9/2]. Here 
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is the action of the bare J J and Sg is the generating func- 
tional of current fluctuations of the conductor introduced 
above. In Eq. J3J) (fi r = ti/2e denotes the reduced flux 
quantum, Cj is the capacitance of the JJ, and the tilted 
washboard potential U (9) = —Ej[cos(9) — s 9], where Ej 
is the Josephson energy and s — h/I c - The factor of 2 in 
the argument of Sg arises from the fact that the voltage 
across the conductor equals the voltage Vj = (Ti/e)(9/2) 
across the JJ. 

In the MQT regime the partition function of the iso- 
lated JJ is dominated by the so-called bounce trajectory, 
an extremal SS[6] = periodic path in the inverted bar- 
rier potential. By approximating a well-barrier segment 
of U(9) around a well minimum 9 m by a harmonic+cubic 
potential, V(59) = {MQ 2 /2)S9 2 (1 - 59/S9 ) with 59 = 
9 — 9 m , one finds an analytic solution in the limit of van- 
ishing temperature, i.e., 
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Here, ft = ft(s) is the frequency for small oscillations 
around the well bottom (plasma frequency) at bias cur- 
rent s, M = (filC j, and <5#o denotes the exit point deter- 
mined from U(9m) = U(9 m + 89q). The corresponding 
MQT rate reads 



36 V b 



T = 6y/6SlV b /hn exp ( -y^ 



where V b = (2A/fi 2 /27) <56>g is the barrier height . 

Following the theory of the effect of an electromagnetic 
environment on MQT |12| , the partition function can now 
be calculated for arbitrary coupling between detector and 
conductor based on a numerical scheme developed in . 
Analytical progress is made when the noise generating 
element has a dimensionless conductance g -C Ej /Ml so 
that the influence of the noise on the MQT rate can be 
calculated by expanding about the unperturbed bounce 
which gives 
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The correction Sg[9b /2] is usually dominated by the sec- 
ond cumulant C2 (width) and the fourth cumulant C4 
(sharpness). Note that this approximation still contains 
the full dynamics of detector and noise source since any 
approximation relying on a time scale separation, as e.g. 
the adiabatic limit considered in |9( , is usually not appli- 
cable. 

Now, in case of a tunnel junction as noise element 
one finds for S g [9b/2} = S T [9 B /2] from |TJ and @ 
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with 
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By expanding the first exponential and performing the 
Fourier transform for each power of Ob(t) separately, the 
relevant part p{u>) — p(c^>) — 2ttS(lu) reads 
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This way, Eq. (JSJ can be cast into 
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with the coefficients 



Here. 
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with Afcfc' = Ak'k so that 7/.^/ ^ only for k + k' 
even. This means that all odd cumulants of the fluctuat- 
ing current vanish according to a vanishing net current 
(I(t)) = through the conductor. Specifically, one finds 
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The terms in the sum (JJJ related to a contribution of or- 
der <5(9g +fc determine the impact of the (fc+fc')th-moment 
of the current fluctuations of the tunnel junction onto 
the MQT process. Since in Eq. JJJJ the term of order 56q 
contains contributions centered around uRiO,fl..., kfl, 
the influence of the (fc + fc')th-moment results from mode 
mixing between fluctuations with frequencies I f2 and I' Q 
where I <k,V < k' . 

In lowest order, k + k' = 2, one gains from Eq. (JJJ the 
Gaussian noise contribution providing a correction to the 
bare MQT rate 



1 j. = 1 exp 
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Apparently, this reflects the well-known fact that Gaus- 
sian noise leads to a reduction of the tunneling rate ^3 • 
At order S9q the sum {7J) gives three contributions, 
namely, k = 1, k' = 3 and k = 3, k' = 1 with A13 = A31 
as well as k = 2, k' — 2 with A22. This leads to 
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FIG. 2: Switching probabilities out of the zero voltage state 
("s-curve") of a J J in parallel to a tunnel junction (solid) and 
to an ohmic resistor (dashed) with identical second cumulant 
for various pulse lengths of the bias current, from left to right: 
lms, lfis, 10ns. Dotted lines below each pair of s-curves dis- 
play the differences of the corresponding switching probabil- 
ities. Parameters are ^jEj/Ec = 10, g T = 2, Q(s = 0)=100 
GHz. 



so that the fourth order cumulant of the current noise 
contains both, fluctuations that suppress tunneling (re- 
lated to A22) and fluctuations that increase MQT (re- 
lated to A31). Since 2^31 — A22 > 0, the total impact of 
the fourth moment leads to an enhancement of the MQT 
rate. 

To obtain explicit results, one derives for the har- 
monic+cubic potential V(S9) an expression for the am- 
plitude S9q (s) of the bounce, namely, 



S6 {s) = 3 



Accordingly, the barrier height scales with the dimension- 
less current s as V b (s) = (2J5j/3)( l - s 2 ) 3 /2 /s 2 , and the 
plasma frequency reads O(s) = [y/2E~jE~c / h] (1 — s 2 ) 1 / 4 
with charging energy E g = 2e 2 /C Experimentally, in 
the standard procedure [Tjj to measure the MQT rate, 
a bias current pulse of height lb and duration t is adia- 
batically turned on and a voltage pulse is detected when 
the JJ switches to its finite voltage state. This procedure 
is performed a few thousand times to built up switching 
histograms that determine the switching probabilities 

P(s) = l-e- r(s)t . 

In Fig. [21 these so-called s-curves are shown for various 
values of the duration of the current pulse t. The pa- 
rameters chosen are accessible in realistic experiments. 
Apparently, for short pulses when the switching occurs 
for values of the bias current close to the critical current 
of the JJ, the effect of the non-Gaussian noise fluctua- 
tions is completely suppressed due to a decreasing am- 
plitude SOo (s) of the bounce. However, for longer pulses 
(or equivalently, for shorter pulses in the tails towards 
lower s values) the fourth order cumulant leads to a sub- 
stantial influence, mostly dominated by a shift to smaller 
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FIG. 3: B(s) = {-ln[r(s)/a]} 2/3 vs. the dimensionless bias 
current s for a tunnel junction (solid) and an ohmic resistor 
(dashed) with identical second cumulant. In the range of 
larger values for s (inset) a straight line (dotted) is fitted to 
the data. This reference line almost coincides with the data 
for Gaussian fluctuations (dashed) while non-Gaussian noise 
(solid) leads to a pronounced reduction of B(s). Parameters 
are the same as in Fig. |2] 

s-values compared to an Ohmic resistor with identical 
Gaussian noise contribution. As seen from Eq. © this is 
due to an effective decrease of the barrier height by the 
non-Gaussian fluctuations of the tunnel junction. The 
contrast in the switching probabilities between MQT in 
presence of purely Gaussian noise and in presence of a 
tunnel junction is larger than 10% for values of the pulse 
height around s = 0.73. This reflects the exponential 
sensitivity of the JJ in the MQT range even to weak 
non-Gaussian fluctuations. 

For the on-chip detection circuit proposed here, it may 
be difficult to see the impact of the fourth order cumu- 
lant directly in the s-curves since reference curves for 
purely Gaussian noise and the same JJ are not avail- 
able. In the standard analysis of MQT data, the rate is 
obtained by exploiting the s-dependence of the bounce 
action. Namely, it turns out that the function 

B( S ) = {-\n{T( S )/a q }y/ 3 



with the prefactor a q — 6^/ WWbjTi/K gives rise to an es- 
sentially straight line, thus reflecting the s-dependence 
of Vfe. This holds true also for weak purely Gaussian 
noise, but any non-Gaussian contributions to the MQT 
rate can be expected to lead to deviations from this 
scaling behavior. For this purpose, s-curves for shorter 
pulses, i.e. larger s-values, where non-Gaussian effects 
are absent, can be used to determine the slope b of 
B(s) - B(s*) = b(s - s*) (reference point s* close to 1). 
For longer pulses corresponding to smaller values of the 
switching current, this straight line b(s — s*) is then com- 
pared with the actual B(s) as depicted in Fig. EI Purely 
Gaussian noise shows basically no deviations from the lin- 
ear behavior over a wide range of s-values, while closer 
inspection of the analytical rate expression (jHJ reveals 
that small deviations lead to a slight increase of the ac- 
tual -B(s) towards lower s. In contrast, non-Gaussian 
noise displays substantial deviations already for much 
larger values of s and leads to a decrease of the actual 
B(s) compared to the reference line extracted for s close 
to 1. Thus, plotting B(s) allows to directly discriminate 
the impact of higher than second order cumulants in the 
noise fluctuations of the conductor. Further, by fitting 
B(s) with Eq. (|HJ), the coefficient 2A 3 i — A 2 2 related to 
d in Eq. © can be extracted. 

The formulation developed above in Eq. (0} is com- 
pletely general and applies to any nanoscale conductor 
in parallel to a J J. A systematic expansion of the action 
<Sg[#b/2] in powers of (6*s/2) n determines the dynami- 
cal impact of the nth order cumulant C„ onto the MQT 
process. 

To summarize we have proposed a nanoelectrical cir- 
cuit where a JJ placed in parallel to an arbitrary conduc- 
tor acts as detector for non-Gaussian current noise. Since 
no net current flows through the noise source, heating ef- 
fects are suppressed and one obtains access to the even 
order cumulants of the distribution function which are 
notoriously difficult to detect. For experimentally realis- 
tic parameters we have explicitly shown how the fourth 
order cumulant of a tunnel junction can be extracted. 
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knowledged. JA is a Heisenberg fellow of the DFG. 



[1] For a recent review see: Y.M. Blanter and M. Biittiker, 

Phys. Rep. 336, 1 (2000). 
[2] L.S. Levitov, H.B. Lee, and G.B. Lesovik, J. Math. Phys. 

37, 4845 (1996). 
[3] Quantum Noise in Mesoscopic Physics, edited by Y.V. 

Nazarov, NATO Science Series in Mathematics, Physics 

and Chemistry (Kluwer, Dordrecht, 2003). 
[4] B. Reulet, J. Senzier, and D.E. Prober, Phys. Rev. Lett. 

91, 196601 (2003). 
[5] R.K. Lindell et al. , Phys. Rev. Lett. 93, 197002 (2004); 

E.B. Sonin, Phys. Rev. B 70, 140506(R) (2004). 
[6] T.T. Heikkila, P. Virtanen, G. Johansson, and F.K. Wil- 

helm, Phys. Rev. Lett. 93, 247005 (2004). 
[7] J. P. Pekola et al, cond-mat/0502446 (unpublished). 



[8] J. Tobiska and Y.V. Nazarov, Phys. Rev. Lett. 93, 

106801 (2004). 
[9] J.P. Pekola, Phys. Rev. Lett. 93, 206601 (2004). 

[10] B. Reulet, cond-mat/0502077 (unpublished). 

[11] G. Schon and A.D. Zaikin, Phys. Rep. 198, 237 (1990). 

[12] A.O. Caldeira and A.J. Leggett, Ann. Phys. (New York) 
149, 374 (1983). 

[13] H. Grabert, P. Olschowski, and U. Weiss, Phys. Rev. B 
36, 1931 (1987). 

[14] A finite capacitance Ct of the tunnel junction can easily 
be taken into account by replacing Cj by C = Cj + Ct- 

[15] M.H. Devoret et al. in: Yu. Kagan and A.J. Leggett 
(eds.), Quantum Tunneling in Condensed Media, (Else- 
vier, Amsterdam, 1992), p. 313. 



